Non-adiabatic transitions in a non-symmetric optical lattice 



O 

C 
3 



43 

Oh- 
i , 

> 

a* 



> 
m 

o 

O 

o 



L. Morales-Molina and S. A. Reyes 
Departamento de Fisica, Facultad de Fisica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 

We study Landau-Zener interband transitions for a non-symmetric optical lattice in the presence 
of an external force. We show that gain and losses of the light beam, as well as the relative occupation 
probabilities of the bands involved in the transitions can be accurately managed upon tuning the 
amplitude of the non-Hermitian component of the lattice. Exact expressions for the transition and 
non-transition probabilities for a non-symmetric system obtained within a two-mode approximation 
are provided. These equations successfully account for the main features of the transitions in the 
optical lattice. We also interpret the non-conventional Bloch oscillations at criticality studied in 
Phys. Rev. Lett. 103, 123601 (2009) as a series of a Landau-Zener transitions. 
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I. INTRODUCTION 

Quantum systems described by non-Hermitian Hamil- 
tonians with real spectra have drawn great interest in re- 
cent years These systems are invariant under parity 
and time reversal symmetry (P7~-syrnmetry) and their 
behavior is quite different from Hermitian ones. Impor- 
tant issues, such as definition of an inner product, uni- 
tarity, as well as the definition of observables for this 
new "quantum mechanics" have been intensively studied 
Q. As a result, new avenues have been opened for the 
observation of exotic phenomena associated with non- 
Hermiticity. 

In particular, recent developments in optics have al- 
lowed the construction of "complex crystals" with unique 
properties such as violation of Friedel's law Q. These 
systems are realized in waveguide arrays that include 
gain and loss regions. Furthermore, a restricted class 
of complex crystals may be built so that they posses the 
aforementioned "PT-symmetry. Along these lines, very 
recent experiments realized for two coupled waveguides 
have revealed a new behavior for light propagation in 
■PT-symmetric optical systems 

Additionally, driven periodic complex potentials have 
proven to be a fertile ground for the study of novel 
phenomena in wave mechanics such as non-conventional 
Bloch oscillations In this paper, we consider the sce- 
nario in which a DC force is applied to a PT-symmetric 
complex lattice. As in conventional quantum mechanics, 
the application of an external force may lead to inter- 
band transitions and interference phenomena. However, 
since in this case the time evolution does not obey the 
unitarity condition, such transitions are fundamentally 
different. Motivated by this idea, in the present work we 
provide a detailed analysis of the Landau-Zener transi- 
tions for a light beam propagating in a complex potential. 

For this purpose simulations of the interband transi- 
tions are performed revealing that the behavior of the 
relative occupation probabilities for the bands can be ac- 
curately managed upon tuning the amplitude of the non- 
Hermitian component of the lattice. In particular, the 
power can be increased or decreased during the transi- 



tion process depending on the sign of the non-Hermitian 
part of the potential. 

Complementary to our numerical results, we develop 
an analytical framework that describes the transitions 
near the level crossings within a two-mode aproxima- 
tion. Notably, the analytical results not only provide 
deep insight into the dynamics of the system, but also 
accurately reproduce the main features observed in our 
numerical simulations. As a special case, we interpret the 
non-conventional Bloch oscillations studied at criticality 
in [J] as a sequence of Landau-Zener transitions. 

The paper is organized as follows: In section II, we de- 
scribe the model for light propagation in periodic struc- 
tures with complex refractive index within the paraxial 
regime. Next, in section III we consider a two-mode ap- 
proximation for the description of the dynamics in the 
complex lattice and provide analytical expressions for its 
non-adiabatic behavior. In section IV, we present a de- 
tailed explanation of our numerical results for the sys- 
tem considered in section II and compare them with the 
theoretical predictions obtained within the two-mode ap- 
proximation. In section V we summarize our results and 
discuss prospects for future developments. 



II. THE MODEL 

Within the paraxial approximation the propagation of 
a monochromatic light beam with wavelength A through 
the lattice can be described by the Schrddinger type equa- 
tion @ 



iXdzty 



X 2 d 2 ^ 
2n s dx 2 



+ U(x)ip - Fx = n tp - Fxip (1) 



where X = X/2n, n s is the refractive index of the sub- 
strate. The potential is given by 



U(x) 



n{x), 



(2) 



where n(x) — n(x + a) is the periodic refractive index 
profile and a is the periodicity of the lattice. The time- 
like coordinate Z is the propagation distance along the 
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waveguides and x is the corresponding transverse axis. 
In the last term of ((T|), F represents the gradient of the 
refraction index that mimics the effect of an external DC 
force. 

For a complex potential, the imaginary part of U (x) 
renders H into a non-Hcrmitian Hamiltonian. In the 
following we shall consider U(x) = U\Cos(2ttx / a) + 
iU 2 sin(27ra;/a), which satisfies U(x) — U(—x)* and 
thus Ho is PT-symmetric 4]. The spectrum of Ho 
has real values for U 2 < f/| rjt = U\, otherwise the 
bands merge giving place to pairs of complex conjugate 
values Q. Henceforth, we consider U 2 < U 2 nt , un- 
less otherwise specified. The eigenfunctions of Ho are 
the Floquct-Bloch band mode functions <f>q, n (%) with 
Bloch-momentum q which fulfill the secular equation 
Ho4> q ,n(x) = E n (q)4> q ^ n (x), where E n (q) are the eigenen- 
ergies corresponding to the n-th Floquet-Bloch band. 
While as in real crystals E n (—q) = E n (q), the wave- 
functions under the parity transformation x — > —x and 
q — > —q obey the general relation <f)\ n (x) — 0_ gjn (— x), 

where 4> qn are eigenfunctions of H . Moreover, the or- 
thogonality condition becomes ( l>q,n( x ) ( i ) q' ,m(x)dx = 
d q , n 5{q - q')S nim with d q , n = ±1 @. 

It is important to notice that for the 'PT-symmetric 
Hamiltonian Ho the power of the propagating light beam 
remains constant provided U 2 < U£ nt . If in addition we 
consider the presence of an external force as in Eq. (JTJ , 
the 'PT-symmetry is broken and the power is not longer 
a conserved quantity. Thus, under the action of the force, 
transitions to other bands may occur which do not con- 
serve the power. Recently, the Gaussian wavepacket dy- 
namics of this system has been analyzed in the presence 
of a small external force Q . In the following sections we 
generalize the analysis for arbitrary values of F. 

To gain understanding of the interband transition pro- 
cesses, we consider as input excitation a plane wave with 
normal incidence. Hence, the wavefunction can be writ- 
ten as a superposition of the Floquet-Bloch eigenfunc- 
tions, i.e. %j) q {x,Z) = E„ c n(Z)</>«(x), where cj> q (x) = 
w q l (x)e lqx and w^(x + a) = w q l (x). Expanding w%(x) in 
a Fourier basis, we can write 

%{x, Z) = Y^ a q (Z)e t{2kl+q)x , with k = n/a 
1 

and a?(Z) = J2<Kn- ( 3 ) 

n 

Here, b\ n with I = 0,±1,±2... are the coefficients of the 
Fourier expansion for w q (x). Substituting this into Eqs. 
(Q312]) we get the dimensionless equation for the dynamics 
of af, 

8a q 

%-£ = (21 + qfaf + (V, + V 2 )af +1 + (Vi - V 2 )af_ x (4) 

where z = Z£k/X, q = q/k and Vj — Uj/2£k, j = 1,2 
with Ek — X 2 k 2 /2n s . It is convenient to write the evolu- 
tion equation in the usual quantum mechanics notation 

id x \a«(z))=H\a«(z)), (5) 
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FIG. 1: (Color online) Energy spectrum f2 = EjEu vs. Bloch- 
momentum q for Vi = 0.2. (a) V 2 = 0.15; (b) V 2 = 0.2. 



where \a q (z)) is the vector with components af and 

H = {H nA = {2l+q) 2 5 n . l +(V 1 +V 2 )5 n . l - 1 + (V 1 -V 2 )S nM1 }. 

(6) 

This Hamiltonian matrix is real but non-symmetric and 
has a real spectrum for values V\ > V 2 . From this expres- 
sion it is also clear that changing the sign of the ampli- 
tude of the imaginary part of the potential V 2 is equiv- 
alent to a transposition of H. Note that the problem 
becomes analogous to the dynamics of a particle moving 
on a discrete chain with non-symmetric hopping ampli- 
tudes (Vi + V 2 and Vi — V 2 ) and on-site energies given by 
(21 + q) 2 . Clearly, for positive values of V 2 , the system 
tends to evolve towards modes with higher energies and 
viceversa. This is a manifestation of the fact that there 
is energy gain (loss) for positive (negative) values U 2 in 
the optical system (TT|). 

A numerical diagonalization of @ results in the energy 
spectrum shown in Fig. [TJ The diagram of energy vs. 
Bloch-momentum reveals that increasing V 2 causes an 
approaching of the bands at the avoided crossings leading 
to degeneracies at V 2 = Vi. 
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III. TWO-MODE SYSTEM 

Let us focus now on the analysis near the avoided 
level crossings that appear at the Bragg scattering points 
q = (2j + 1) with j = 0,±1,±2.... For a small energy 
separation between the modes involved in a given cross- 
ing, one can safely neglect the effects due to any other 
bands. A suitable description can then be provided in 
terms of an effective two-mode Hamiltonian. In particu- 
lar, from the matrix ([6]) and considering the two lowest 
bands it is possible to obtain the corresponding effective 
Hamiltonian 



H 



A+S 
c 2 
A-<5 

2 



— r 



(7) 



where e = 2(1 - 21 - q), A = 2V X and 5 = 2V 2 . 
Unlike the Hermitian case, in Eq. J7| H± t 2 ^ EL^ \ 
and the eigenvalues can therefore have complex val- 
ues. By solving the secular equation for ([7]), we get 

E± = ± \J e 2 + ^ — j — - so that for high values of e the 
spectrum is real. At e — the gap between the energy 
levels becomes E + — E_ = VA 2 — S 2 . On the other 
hand, when S 2 > A 2 purely imaginary eigenvalues arise 
in the region e 2 < (5 2 ~ A 2 )/4, indicating that the system 
is beyond the critical point where the 'PT-symmetry is 
broken. 



A. Landau-Zener for a non-symmetric two-mode 
system 

The presence of a gradient in the refraction index in 
Eq. (JT]) mimicking a DC force causes a linear variation 
of the Bloch-momentum along the propagation distance, 
i.e. q — az, where a is proportional to the force. As a 
result the diagonal elements of matrix (JT)) are modified as 
e = —2az + €q, where eo = 2(1 — 21). By introducing the 
change of variable t = z — eo/2a we obtain e = — 4r where 
f3 = 4a. Thus, at t = the energy levels are only sep- 
arated by the gap. In Landau-Zener theory the variable 
t represents time and (3 is the speed for the variation of 
e whose tuning controls the transitions between the two 
bands [7j, 8| . To study the interband transition, the initial 
state is prepared away from the avoided crossing where 
|/3f| 3> \/A 2 — 8 2 . In the conventional Landau-Zener 
framework this amounts to a preparation at t — > — oo. 
For simplicity we consider |ai(— oo)| 2 = 1 for the lowest 
energy state and 02 (—oo) = for the excited state. 

Now, as the system evolves, we find that after suffi- 
ciently long time the asymptotic values for the intensities 
are given by @ 



|02(00)| S 



P 



exp 



2/3 



for the excited state and 



K(oo)| 2 = r(i-p) 



(8) 



(9) 




FIG. 2: (Color online) Schematic representation of two energy 
levels in the limit \S\ = A. (a) 5 = A; (b) S = -A. 



for the ground state, where F = (A + S)/(A — 6). 

Since for this two mode system the power is given by 
p = 1 0,1 1 2 + |a.2 1 2 , after the transition we have that p = 
r(l— P)+P, which in the limit of very slow variation (3 <C 
(A 2 — S 2 ) reduces to p — \ai\ 2 =T. Thus, permutation of 
the nondiagonal matrix elements in (|7|) results in power 
variations of the order of 4A<5/ (A 2 — 15 2 ). This shows that 
even a slow navigation through the spectrum can greatly 
affect the power due to the coupling with the upper state. 
In this regard, the quantities |o.2 1 2 and |ai| 2 can only be 
referred to as the relative transition and non-transition 
probabilities, respectively. 

To get some insight into Eq. ©, we analyze several 
regimes for the parameter S. First of all, if we set S = 
into the equations above, corresponding to a real poten- 
tial (V2 = 0), the original Landau-Zener probability tran- 
sition Plz — exp[— 7rA 2 /2/3] is recovered and the non- 
transition probability becomes (1 — Plz)- If < <5 < A, 
the relative transition probability P increases for increas- 
ing values of S. The same holds for the quantity T(l — P), 
even though the quantity (1 — P) is decreased. On the 
other hand, for — |A| < 5 < and increasing \S\, the 
relative nontransition probability decreases. 

Consider now the limiting case S — > A with /3 ^ 
where the gap vanishes, i.e. there is a band crossing. In 
this regime the expressions given in reduce to 



lim |oi| 

5->A 



2tt- 



A. 



"2 



1. 



(10) 



Observe that for small values of (3 the intensity in the 
lowest mode is greatly increased for any A > 0. In con- 
trast, for negative <5, the asymptotic intensities become 



lim I ai| 

«->•- A 



0. 



«2 



1. 



(11) 



It is interesting to note that, contrary to the conven- 
tional Landau-Zener transition process for a vanishing 
gap, there is a non-zero occupation probability for the 
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FIG. 3: (Color online) Power (p) as a function of z. \a\ = 
1 x 10 -3 and V 2 = 0.1. solid line:a > 0, dashed line: a < 0. 



ground state, as observed in the first scenario (fTU|) . On 
the other hand, the relative transition probability to the 
excited state is |a2 1 2 = 1 (see Fig. (JS^,)). Remarkably, in 
the second scenario, there is a complete transition to the 
excited state with |ai| 2 = (see Fig. (Jib)) as expected 
for a standard Landau-Zener process. The fundamental 
differences between these two situations are a manifesta- 
tion of the non-symmetric form of the system (J7J . 

So far we have shown the results for a non-adiabatic 
transition for a single two-level system. Based on the 
same scheme, these results can be applied to the descrip- 
tion of a series of avoided crossings in the optical multi- 
level system ^ . Let us now assume that the input exci- 
tation lies completely in the lowest mode, i.e. |aj| 2 = 5±j. 
As the force F is turned on, the bloch momentum is 
modified and the system undergoes a succesion of tran- 
sitions through several (avoided) level crossings. Note 
that due to the lack of coupling among higher bands, the 
only allowed transitions are those between the two lowest 
modes. Consequently, it is straightforward to find that 
for F^O the total power after n-th consecutive avoided 
crossings is given by 

n-l 

/ 9„ = r"(i-p)" + pJ]r(i-p) 1 , (12) 

i=0 

whose expression in the limit 6 — > A reduces to p n — 
^™ =0 A l (cf. Eq. (ITU|) ). Note that the power always 
exceeds unity. On the other hand, taking 6 —> —A the 
power reduces to p n = 1. 

In expression p2[) the first term accounts for the 
change of probability density in the lowest level, while 
the second term represents the part of the probability 
density that is transferred to higher energy levels after n 
avoided crossings. 
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a 

FIG. 4: (Color online) Probability transition P vs. rate a. 
(a) V 2 = 0, (b) V 2 = 0.15, (c) V 2 = 0.19. Dashed lines 
superimposed correspond to the results obtained with Eq. ([5|. 



IV. PT-SYMMETRIC OPTICAL WAVEGUIDE 
ARRAYS 

The optical realization of a medium with VT- 
symmetric properties has paved the way for the synthesis 
of artifical optical periodic lattices with 'PT-symmctry. 
This makes the observation of interband transitions in 
a PT-symmetric lattice a target within reach in a near 
future. 

As explained above interband transitions are generated 
by the presence of a gradient in the refractive index. In 
order to understand these proceeses we look at the transi- 
tion points, i.e. the avoided energy crossings of the band 
structure of the complex lattice. Beam splitting occurs 
at these scattering points due to non-adiabatic transi- 
tions and the intensity of the beam is modified in the 
process. We focus our analysis on the parameter window 
< V2 < V" 2 crlt , where there are no singular or excep- 
tional points. As previously noted, the action of the fic- 
titious force in Eq. ([T]) can be effectively replaced in Eq. 
((3]) by a Bloch-momentum q that depends linearly on z, 
i.e. q = az with a = F/(k£k)- 

As the initial state we choose an input excitation 
that populates the lowest band with Bloch-momentum 
9 = <7im = 0, which can be realized by shining a normal 
incident plane wave in the z direction. We start consid- 
ering first the limit of slow variation, \a\ <C (Vf 2 — V%)- 
Fig. ^p) shows the numerical results for the evolution 
of the power p in this limit for both positive and neg- 
ative a. Interestingly, for a > the asymptotic value 
for the power is p + w 3 whereas for a < is p~ m 1/3. 
Moreover, these values can also be accurately estimated 
using Eq. (112[) in the corresponding limit for a single 
avoided crossing [n = 1), giving p + = 3 and p~ = 1/3 
respectively. 

Next, for a complete analysis of the interband tran- 
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sition process, we consider a range of values of a going 
from the adiabatic limit considered before, to a strongly 
non-adiabatic regime where the energy associated to the 
force becomes much larger than the interband gap. To do 
so, for each a we perform simulations for the evolution of 
the system from = to <jfi n = 1.8 and measure the 
resulting occupation of the excited state. The latter is 
obtained by projecting the evolved wavefunction $ onto 
the eigenstate <$ q 2 corresponding to the second lowest 

band of V) with q = q~f in leading to P = |(<^ i2 |$)| 2 - I n 
FigH] we plot the numerical calculations for the transi- 
tion probability P vs. a. It is observed that the profile 
for P(a) in the non-symmetric case remains very simi- 
lar to the one for V<i = if the system is far from the 
critical point (V2 <C V 2 cnt )- In fact, the difference be- 
tween the non-symmetric and symmetric cases becomes 
noticeable only when V2 is fairly close to its critical value. 
This happens because as V2 approaches V^ 7 " 1 *, the gap is 
drastically reduced, leading to a strong increase in the 
transition probability. As the limit V2 — V± is reached, 
the gap vanishes and the transition probability becomes 
unity for any value of a =/= 0. Using Eq. (j5J) derived 
within the two-mode model, we obtain the correspond- 
ing results which are also plotted in Fig. (QJ. Notably, 
the approximated analytical results match the numerical 
simulations impressively well, complementing the agree- 
ment found in the adiabatic limit. 

As the system is driven through Bloch-momentum 
space, a series of avoided crossings in the spectrum may 
lead to multiple transitions. In this regard, we carried 
out simulations going from qi n i — to q~fi n — 3.9 with 
a = 3 x 10~ 2 , such that the window in the q axis in- 
cludes two avoided crossings. Fig. [5] shows the power p 
as a function of the propagation distance z. The funda- 
mental feature in both panels is the emergence of steps at 
the Bragg scattering points as the power evolves. While 
for a > (Fig. [S^,) the power exhibits growing steps 
after every avoided crossing, in the case of a < the 
power decreases in a sequence of small steps (Fig. [5p). 
To understand these two different scenarios we follow the 
evolution of the beam power through the band struture 
with the help of expression (TT2)) . 

Let us first analyze the critical situation V2 = V\ where 
the two lowest bands cross: For a positive force our an- 
alytical results predict a behavior similar to Bloch oscil- 
lations, reflected in the large value of the intensity of the 
ground state after n consecutive crossings (A n ) and with 
the remaining power (^™Zq A 1 ) transferred to higher en- 
ergy bands ^j. This process is schematically depicted in 
Fig. [TJd by arrows pointing to the right. In contrast, for a 
negative force, the ground state has a null intensity after 
n succesive crossings which clearly indicates the absence 
of Bloch oscillations. Here, a full transition to the ex- 
cited state takes place after undergoing the first crossing 
and since there is zero coupling between higher bands, 
the state follows a parabolic trajectory as schematically 
depicted in Fig. [T}d with arrows pointing to the left. 

For nonzero gap, transitions become determined by the 
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FIG. 5: (Color online) Power p as a function of the prop- 
agation distance z. (a) a > 0; (b) a < 0; In both panels 
solid lines represent: Critical case V2 = V\ (thin red line) and 
V2 = 0.15 (thick black line). Dashed lines mark the asymp- 
totic values calculated with the corresponding Eq. (|12|) . The 
other parameters are Vi = 0.2, |a| = 3 X 10 -2 and the critical 
case is computed with V% = 0.2 — 10 -7 . 



tradeoff between the applied force and the gap, resulting 
in a rich picture of Bloch-oscillation-like behavior com- 
bined with interband transitions whose average power 
can be estimated using expression (jT2J) (see horizontal 
dashed lines in Fig. [5]). Notice that for a > the height 
of the steps is decreased when the amplitude of the imag- 
inary component of the potential is reduced (Fig. EK), as 
expected due to the associated absorption process. Con- 
versely, for a < there is an attenuation of the power at 
each transition as observed in Fig. [SJd. Once again the 
applicability of the theoretical equations derived within 
the two-mode formalism is demonstrated by the excellent 
agreement with numerical simulations. 

We would also like to point out that all of the results 
mentioned above for the two-mode approximation cor- 
respond to asymptotic values far from the avoided level 
crossing. Nevertheless, close to the transition region, the 
exact two-mode calculation exhibits an oscillating behav- 
ior similar to the one present in the simulations shown 
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in Fig. [5j Clearly, these oscillations arise as a result of 
the interaction between the bands and are damped out as 
the coupling becomes negligible when the system evolves 
away from the (avoided) level crossing. 

V. CONCLUDING REMARKS 

In summary we have investigated the Landau-Zener 
transitions for a "PT-symmetric complex lattice that can 
be realized in optical waveguide arrays. Firstly, the anal- 
ysis was performed using an exact calculation for the two- 
mode approximation and later extended to the complete 
multi-mode system driven through a series of avoided 
crossings. Overall, we found that the two-mode frame- 
work provides a simple and transparent understanding of 
the main features of the dynamics in the complex crystal. 

The analytical expressions presented above not only 
reflect the non-reciprocity of the lattice, but also accu- 
rately predict the simulation results including the more 
general scenario with multiple crossings. Also, it is found 
that the power is very susceptible to the system speed 



as it moves accross the Bragg scattering points, where 
the passage time determines the amount of energy ab- 
sorbed/emitted by the light beam. This mechanism can 
be used as a way to control the intensity of the light beam 
in a complex waveguide array. 

In a different context, complex optical lattices for cold 
atoms have been built [10|- In this regard, investigation 
of transitions for matter waves between bands realized 
for real optical lattices flll - [T3l | could in principle be ex- 
tended to PT-symmetric systems. Furthermore, the re- 
sults exposed on this work could eventually serve as a 
benchmark for future developments in the several areas 
that involve Landau-Zener processes and interferometry 
for non-symmetric Hamiltonians. 
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